Absence of Dirac Electrons in Silicene on Ag (111) Surfaces 
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We report first-principles calculations that clarify stability and electronic structures of silicene 
on Ag(lll) surfaces. We find that several stable structures exist for silicene/Ag(lll), exhibiting a 
variety of images of scanning tunneling microscopy. We also find that Dirac electrons are absent near 
Fermi energy in all the stable structures due to buckling of the Si monolayer and mixing between Si 
and Ag orbitals. We instead propose that either BN substrate or hydrogen processing of Si surface 
is a good candidate to preserve Dirac electrons in silicene. 
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Graphene is a hexagonally bonded carbon sheet in 
which sp 2 hybridized electrons (a electrons) form a hon- 
eycomb structure and the remaining tt electrons follow 
the massless Dirac (Weyl) equation. The energy bands 
show the linear dispersion (Dirac cone) at the Fermi level 
(E-p) at particular symmetry points, K and K f , in Bril- 
loine zone (BZ) |l|. This causes fascinating properties 
such as the anomalous quantum Hall effect [2|, l3| and the 
unexpected magnetic ordering 0, 0]. Moreover, the in- 
plane high mobility intrinsic to graphite and its robust- 
ness place graphene as an emerging material to develop 
post-scaling semiconductor technology challenged by the 
cutting edge miniaturization [6]. 

In order to maximize technology connectivity, however, 
a hexagonal sheet consisting of Si which is located just 
below C in the periodic table has obvious superiority. 
Spin-orbit interaction stronger in Si than in C, possi- 
bly inducing the quantum spin-Hall phase [7] , is another 
intriguing factor in science viewpoints. Several experi- 
mental efforts to synthesize such a Si sheet (silicene) has 
been indeed done recently J8[ [s| , and buckled silicene on 
and on ZrB2 [13| surfaces are identified by 
scanning tunneling microscopy (STM) measurements. 

Calculations [14| in the local density approximation 
(LDA) in the density-functional theory (DFT) have clar- 
ified that a hexagonally bonded flat Si sheet is unstable 
to the buckling and the resultant, still hypothetical, free- 
standing silicene with the buckling of 0.44 A shows the 
Dirac cone as in flat graphene. It is therefore extremely 
important to reveal whether Dirac electrons exist in the 
real silicene on the substrates. 

Experimentally, the Si layer on the Ag(lll) surface 
has been reported to exhibit several structures with the 
super-periodicityof 4x4 \/l3 x \/l3 [111, [Hj, 

2\/3 x 2V3 [H GJ, [l6| and y/3 x >/3 with respect to 
1 x 1 of the Ag(lll) surface. The former two have been 
well identified, whereas the latter two are still of contro- 
versy in terms of lattice mismatching [10|, • Angular- 
resolvedphotoelectron spectroscopy (ARPES) measure- 
ments [10| suggest the existence of the linear-dispersion 
band near E^ but the spectra are not conclusive. 

In this Letter, we report extensive calculations for sil- 



icene on Ag(lll) surfaces with the periodicity of 4x4 
and Vl3 x \^13 as well as 2\^3 x 2\^3. We find six dis- 
tinctive stable structures with the total-energy difference 
ranging in 69 meV per Si atom. Calculated scanning 
tunneling microscopy (STM) images of the most stable 
structures are in excellent agreement with the experi- 
ments available. An important finding common to all 
the (met a) stable structures is the interaction between sil- 
icene and Ag substrate, leading to the disappearance of 
the Dirac electrons near Ep. We instead propose that 
the hexagonal BN (h-BN) substrate and the hydrogen- 
processed Si [18] are promising to preserve Dirac elec- 
trons in silicene. 

Calculations are performed by our real space DFT (RS- 
DFT) code [19]. The LDA is used to the exchange- 
correlation energy, and the norm-conserving pseudopo- 
tentials are adopted [20j to describe the electron-ion in- 
teraction. The grid spacing in the real-space calculations 
is taken to be 0.16 A corresponding to a cutoff energy 
of 108 Ry in the plane-wave-basis-set calculations. The 
silicene-covered Ag surface is simulated by a repeating 
slab model in which five- atomic- layer slabs are separated 
from each other by a 14- A vacuum region. The geome- 
try optimization is performed until the remaining forces 
become less than 26 me V/A dH. 

We have first performed extensive search for stable ge- 
ometries of silicene on Ag (111) surface. It is recognized 
that the 3x3 super-periodicity of silicene is commen- 
surate to the 4 x 4 Ag(lll) with the lateral Si spacing 
of 2.23 A, whereas the x y/7 periodicity of silicene 
is to the vT3 x ^13 and the 2V^ x 2V3 Ag(lll) with 
the lateral Si spacing of 2.28 A and 2.19 A, respectively. 
We have found two distinctive structures each for the 
4x4, the Vl3 x Vl3 and the 2\f?> x 2\f?> periodicities. 
The structural parameters of obtained stable (S) and 
metastable (MS) structures are shown in Table U along 
with the corresponding adsorption energy i? a d, which is 
denned by E ad = (E Ag + A^si^Si - E Si/Ag )/N S i. Here 
i^Si/Ag and E Ag are the total energies of silicene/ Ag(l 11) 
and clean Ag(lll), respectively, and Nsi is the number 
of Si atoms in the silicene. As for the chemical potential 
/isi of Si, we use the total energy of an isolated Si atom. 
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TABLE I. Calculated adsorption energy E a d (eV/Si-atom) 
and structural parameters of the stable(S) and the metastable 
(MS) structures for silicene on Ag(lll) surfaces with the pe- 
riodicities of 4x4, ^13 x Vl3 and 2^3 x 2^3. The cfei-Si and 
^Ag-si are the average distances (A) between Si atoms and 
between bottom Si and Ag surfaces. The Az ( A) shows the 
Si coordinates in the surface-normal direction relative to that 
of the bottom Si. 



Periodicity 


str E ad 


dsi-si 


0?Ag-Si 


Az 


4x4 


S 5.972 


2.29 


2.13 


— — 0.78 




MS 5.968 


2.29 


2.11 


0.57, 0.70, 0.90 


Vl3 x Vl3 


S 5.975 


2.32 


2.13 


0.36, 0.59, 0.79 




MS 5.971 


2.31 


2.13 


0.19, 0.52, 0.90 


2a/3 x 2^3 


S 5.922 


2.27 


2.16 


— , 0.36, 1.19 




MS 5.906 


2.27 


2.21 


— , — , 1.06 



Calculated adsorption energies for the S and MS struc- 
tures for the Vl3 x \/l3 and the 4x4 surfaces are close 
to each other: E a( ± is largest for the S structure of the 
\fl3 x vT3, followed by those for other three structures 
with less than 10-meV decrease. It is of note that E^ of 
the silicene is comparable with the cohesive energy of Si 
in diamond structure (6.02 eV), and larger than that of 
free-standing silicene (5.27 eV in our calculations). The 
adsorption energies for the 2^/3 x 2^/3 surface is relatively 
smaller by 50 - 70 meV, reflecting the compression of the 
Si layer on this surface. The energetics obtained is con- 
sistent with that the 2\/3 x 2\^3 surface is less frequently 
observed [ToHil, E3, E3 • 

One of the characteristics common to all the S and MS 
structures is the buckling of silicene: The surface-normal 
coordinate z of Si has 2-4 different values depending on 
the structures, as shown in Table |U The amount of the 
buckling is larger on average than that of free-standing 
silicene. The average of the Si-Si distance dsi-Si is larger 
than that of free-standing silicene (2.24 A in our calcu- 
lation). It is of note that the averaged distance between 
Si and Ag surface infers the formation of chemical bonds 
between silicene and Ag. 

The complex pattern of the buckling we find leads to 
a rich variety of STM images. Fig. [1] shows calculated 
STM images [22] of the S and MS structures for the sil- 
icene/ Ag with 4 x 4, V^xV^, and 2 V3x 2^/3 periodic- 
ities. The calculated STM image for the S structure with 
each periodicity is in excellent agreement with the exper- 
iments EMI 11611. Our S structure for the 4x4 surface 
seems essentially identical to the structures proposed in 
loL 11 1 . The structures with other 
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periodicities have been proposed 
ment with the experimental STM images is poor. On the 
contrary, the calculated STM images in the present work 
reproduce characteristic features in the experiments: The 
moderately bright region of the lower half of the lateral 
cell in the y/l3x V^13 structure [Fig. ffl^c)] and the super- 
hexagonal structure in the 2 y/3 x 2\^3 structure [Fig. 




FIG. 1. (color online) Calculated STM images (right parts) 
and lateral atomic positions of Si atoms (left parts) of stable 
structures for silicene/Ag(lll) surfaces. The S (a) and the 
MS (b) structures for the 4x4 surface, the S (c) and the 
MS (d) structures for the y/l3 x y/l3, and the S (e) and the 
MS (f) structures for the 2^3 x 2^3. The buckling of the 
silicene is represented by the size of each ball in the left parts: 
The largest red and smallest blue balls depict the top and the 
bottom Si atoms and other balls do the Si atoms between. 
The lateral unit cells are indicated by the orange lines. 



Ute)]. The calculated adsorption energy of each struc- 
ture corroborates the identification proposed here. The 
MS structure with each periodicity exhibits the STM im- 
age different from that of the S structure [Figs. H(b), (d), 
(f)]. These images are not observed so far, but are likely 
to be observed in future since the difference in be- 
tween the S and the MS structures is small. 

We are now in a position to discuss the electronic struc- 
ture. Figure [2] shows the energy bands of the S structures 
of the 4x4 and the \/V3x a/13 silicene/ Ag(l 11) surfaces. 
We also show the energy bands of the isolated 4x4 and 
the a/13 x V^13 silicene layers which are peeled from the 
Ag substrate in the S structures. In the low-buckled free- 
standing silicene, the energy bands clearly show Dirac 
cone and the corresponding Kohn-Sham (KS) orbitals at 
K point have characters of the bonding tt and the anti- 
bonding 7r (labeled as 7r*), as is shown in Fig. [3ja). In the 
isolated silicene peeled from the Ag substrate, the tt and 
7r* states still exist near Ep [Figs. [2fc) and (d)]. They 
are almost degenerate at Ep for the x \fY3 struc- 
ture, whereas they split for the 4x4 structure, opening 
a gap of 0.3 eV. It is recognized that the high buckling 
opens a gap due to incorporation of sp 3 hybridization 
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FIG. 2. (color online) Calculated energy bands of the S struc- 
tures for the 4 x 4 (a) and the Vl3 x Vl3 (b) silicene/Ag(lll) 
surfaces in the lateral BZ. The K point in the BZ of 1 x 1 sil- 
icene is folded on T in the 4x4 and on K in the \/l3 x \^13 sil- 
icene/Ag(lll) surfaces, respectively. Ef is shown by dashed 
lines. Calculated energy bands of the isolated silicene layer 
peeled from the S structure of the 4 x 4 (c) and the vT3 x a/13 
(d) surfaces are also shown. States marked by squares have 
characters of tt and tt* , or the mixed tt and tt* (see text). The 
linear band of mixed tt (it*) state in (a) is indicated by the 
black dots. 



FIG. 3. (color online) Contour plots of KS orbitals of the 
mixed tt (tt* ) states in silicene/Ag(lll) on the plane perpen- 
dicular to the silicene sheet, (a) KS orbitals (absolute value) 
of the tt (upper panel) and tt* (lower panel) states for the 
freestanding silicene at K point in lxl BZ. (b) KS orbitals 
of the four mixed tt (tt*) states marked in Fig. [2ja) in the 
4x4 silicene/Ag(lll) at F point, (c) KS orbitals (absolute 
value) of the two mixed tt (tt*) states marked in Fig. [2jb>) in 
the vT3 x vT3 silicene/Ag(lll) at K point. The blue and 
white balls depict Si and Ag atoms, respectively. The contour 
values are normalized to the maximum. 



13l Il4j . The energy bands of the isolated 4x4 silicene 
is the case (see Az in Table [I]). In the isolated y/l3 x y/l3 
silicene, the gap is very small presumably due to uneven 
arrangement of Si atoms with small Az. 

In real silicene/ Ag(l 11), the tt and tt* states disappear 
near Ep. Detailed analyses of the Kohn-Sham orbitals 
clarify that the states which have the components of tt 
or tt* orbitals exist in the valence bands. They are at 
-1.11 eV (-1.28 eV) below E F for the 4 x 4 (the Vl3 x 
y/l3) structure [Figs. [2] (a) and (b)]. Our calculation 
shows that the amount of the electron transfer from Ag 
surface to silicene is less than 0.43 (0.36) electron per 
unit cell for the 4x4 (a/I3 x \/l3) surface. This electron 
transfer causes the shift of Ep to the higher energy bands 
of silicene by about a few tenths of eV, which is much 
smaller than the calculated shift of possible Dirac points. 

Figures [3fb) and (c) show KS orbitals of the above 
states. It is found that those tt and tt* states strongly 
mixed with Ag orbitals at the top two layers [mixed tt 
and tt* states] . This mixing leads to the formation of Si- 
Ag bonds in the S structures, rendering the states deep 
in the valence bands. The anti-bonding counter parts 
indeed exist around 1 eV above Ep with the bonding- 
antibonding splitting being more than 2 eV. Calculated 
total energy gain, - (E Si/Ag -E Ag -^silicene), is about 0.70 
eV/Si-atom in both the 4x4 and Vl3 x \^13 surfaces. 
Dirac states are substantially modified and shifted from 
Ep due to the mixing of the Ag substrate. We have 



confirmed that the absence of Dirac electrons found here 
is common to all the (met a) stable structures of silicene 
on Ag(lll). 

In order to clarify the reason of the absence of Dirac 
electrons, we peel the silicene layer from the 4x4 Ag 
substrate, as in Fig. E](c), and then put hydrogen atoms 
below the silicene with a fixed Si-H distance of 1.50 A 
as substitutes for the Ag atoms. The energy bands and 
the KS orbitals of the resultant 3x3 silicene/Hig are 
shown in Figs. |4] (a) and (c). We find that the mixed 
tt (tt*) states appear around -3.1 eV at T but the degen- 
eracy is lifted by 57 meV similar to that in Fig. [2] (a). 
The bonding between Si and H atoms renders the mixed 
tt (tt*) states far from Ef- Energy bands of the silicene 
peeled from \^13 x V^13 Ag substrate with the H atoms 
attached show the essentially identical feature. We next 
make the H-attached 3x3 silicene/Hi8 flat and examine 
the band structures. The gap between the lifted mixed tt 
(tt*) states located at around -3.1 eV becomes 25 meV in 
this case, recovering almost the linear dispersion. This 
small gap of 25 meV is due to the lateral distortion of 
Si hexagons in silicene. We finally remove H atoms from 
the above flat silicene 3x3 silicene/Hig. It is then found 
that the mixed tt and tt* states lose the character of hy- 
drogen orbitals and shift upward at Ef [FigHfb)], with 
Dirac electrons recovered. 

We now recognize that two factors are important in 
disappearance of Dirac electrons in silicene: the buckling 
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FIG. 4. (color online) Calculated energy bands of hydrogen 
attached silicene 3 x 3/His (a) and of flattened silicene 3 x 
3 without H (b). The Kohn-Sham orbitals for the four mixed 
tt and tt* states marked by the squares at Y in (a) are shown 
on the plane perpendicular to the silicene sheet in (c). The 
Dirac point comes up at Ef in (b) depicted as a square. The 
blue and pink balls depict Si and H atoms, respectively. 



of the silicene layer and the mixing with the substrate or- 
bitals. Our calculations show that the buckling induces 
the mixing of the tt and the a orbitals and usually breaks 
the hexagonal symmetry, thus breaking the linear disper- 
sion and the lifting of the degeneracy. When the silicene 
stands alone, the tt and tt* states, whichever is mixed or 
not with the a, stay near Ef- Then when the silicene 
is placed on the substrate, the mixing with the substrate 
orbitals induces the mixed states which are shifted down- 
wards, being located deep in the valence bands. Inter- 
estingly, the mixing between tt (tt*) and the substrate 
orbitals have another effect: In silicene 4x4 /Ag (111), 
the mixing between Si and Ag orbitals reduces the mix- 
ing of the 7r and the a orbitals in silicene and recovers the 
linear dispersion presumably due to the partial recovery 
of the hexagonal symmetry [Fig. [I] (a) and Fig. [2] (a)]. 

Finally, we propose other two systems in which the 
character of Dirac electrons of silicene is satisfactorily 
preserved: i.e., the h-BN substrate and the hydrogen- 
processed Si. Figure [5] shows calculated energy bands and 
KS orbitals of silicene on h-BN sheet. In our calculation, 
\/3 x y/3 periodicity of silicene is commensurate with 
\fl x y/7 periodicity of h-BN with the lattice mismatch 
of 0.64 %, and the binding energy between silicene and 
h-BN is 56 meV/Si-atom with the optimized inter-layer 
distance of 3.23 A. The mixed tt and tt* states appear 
near Ef with the tiny gap of 10 meV, exhibiting the 
linear dispersion peculiar to Dirac electrons. 

Figure [6] shows calculated energy bands and KS or- 
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FIG. 5. (color online) Calculated energy bands and contour 
plots of the mixed tt and tt* states at T point of silicene on 
h-BN. Extended energy bands near Ef — are shown in the 
inset. The mixed tt and tt* marked by the square are shown 
on the plane perpendicular to the silicene layer. The blue 
balls depict Si atoms. 




FIG. 6. (color online) Calculated energy bands (a) of 
hydrogen- inserted Si (111) surface [inset of (a)] and contour 
plots of the mixed tt and tt* states (b). The blue and pink 
balls in the inset depict Si and H atoms, respectively. Ex- 
tended energy bands near Ef — are shown in another inset 
and the distributions of the tt (tt*) states at K point marked 
by the square are shown on the plane perpendicular to the 
surface. 



bitals of Si (111) surface with the hydrogen monolayer 
intercalated in the subsurface. The H insertion in the 
subsurface is used to cut Si crystal [18| and the single- 
layer H insertion may be possible in future. We have 
indeed found that such intercalated structure is stable 
with the separation between the top Si layer and the H 
layer of 2.44 A. The binding energy of the top Si layer 
to the H-covered Si (111) surface is 140 meV/Si-atom. 
The top Si layer keeps the hexagonal structure with the 
difference in the surface- normal coordinates Az = 0.48 
A of the Si atoms. Calculated energy bands and the KS 
orbitals clearly show the existence of Dirac electrons at 
E F . 

In conclusion, we have performed first-principles total- 
energy electronic-structure calculations for silicene on Ag 
(111) surfaces. We have found that buckling of the sil- 
icene layer and the mixing of tt and tt* orbitals of silicene 
with the Ag orbitals are prominent. This makes the Dirac 
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electrons of the hypothetical free-standing silicene disap- 
pear near the Fermi level. We have instead proposed that 
either BN substrate or the hydrogen intercalation in Si 
(111) subsurface provide a new stage to preserve Dirac 
electrons in real silicene. 
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